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LINEAR THEORY OF IMPULSIVE VELOCITY CORRECTIONS 


FOR SPACE MISSION GUIDANCE 

By Luigi S. Cicolani 
Ames Research Center 


SUMMARY 


One aspect of midcourse guidance is the determination of the guidance 
correction, usually specified by an equation that gives the required velocity 
change as a function of the current state of the vehicle. With the assumptions 
that impulsive velocity corrections are made and that the equations of motion 
can be linearized about a reference orbit, the required midcourse velocity 
correction will be a linear function of the deviation of the vehicle state 
from the reference state . 

Such linear guidance equations form a particularly simple class which is 
analyzed in this paper from the viewpoints of linear vector spaces and the 
theory of linear equations. Existence conditions and two formulations of the 
guidance equation from the mission constraints are given. Properties of form 
and time -invariant properties of the guidance equations and relations among 
the basic notions of the theory are also discussed. Finally, a number of 
examples from the recent literature are examined. 


INTRODUCTION 


A mission must satisfy certain objectives which, in turn, specify con- 
straints on the vehicle’s orbit. The launch system attempts to place the 
vehicle on a reference orbit that satisfies the constraints, but misses this 
objective slightly. For lunar and interplanetary missions the injection errors 
will usually be unacceptable for mission success and a midcourse navigation 
and guidance system and one or more midcourse corrections will be required. 

The navigation system estimates the vehicle state from observations of stars, 
planets, or the vehicle, and the vehicle guidance system then attempts to cor- 
rect to a satisfactory orbit (one which meets the mission constraints on 
orbital motion) by means of a change of velocity. Because the actual state is 
imperfectly known and the correction imperfectly made, the vehicle is not yet 
on a satisfactory orbit and further observations and corrections may be 
necessary. 

Only one element of the midcourse guidance and navigation process will be 
ex ami ned here - the equation which gives the required velocity correction as a 
function of current vehicle' state . 

In the study of midcourse guidance for the lunar and interplanetary 
missions, it has frequently been adequate to assume that velocity corrections 



are impulsive and departures from a reference orbit are sufficiently small to 
allow the equations of motion to be linearized about the reference orbit. The 
guidance equations that result from these assumptions give the velocity correc- 
tion as a linear function of the deviation of the current vehicle state from 
the reference state. Guidance equations of this type will be termed linear 
impulsive guidance laws herein. A number of such guidance laws designed to 
satisfy various mission constraints have been reported (e.g., refs. 1-5)* All 
of these examples of guidance laws show some similarity of form, suggesting 
that they are mathematically related. It is the object of this paper to deter- 
mine the general properties and logical structure of the class of linear 
impulsive guidance laws. 

The analysis adopts two approaches; in the first we consider the state 
space. The midcourse velocity correction is an attempt to correct some arbi- 
trary state to a satisfactory state, that is, to a state which defines an orbit 
that satisfies the mission constraints. Linear impulsive guidance laws occur 
when the satisfactory states are, collectively, certain types of subspaces of 
the state space, and the guidance law can be formulated from any basis of the 
satisfactory states. Further, such guidance laws have properties of form which 
depend on the dimension of the set of satisfactory states and are invariant 
with time. The set of satisfactory states is, in turn, defined by the mission 
constraints on orbital motion, and a basis of the satisfactory states can be 
generated by a sufficient number of derivatives of the state, evaluated on the 
reference orbits, with respect to appropriate orbital parameters. 

The second approach utilizes the mission constraints directly; from the 
linearized relation between the constraint parameters and the current state, 
the general solution for the velocity correction which satisfies the con- 
straints can be obtained. This construction gives the guidance law from deriv- 
atives of the constraint parameters with respect to the Cartesian components 
of the current vehicle state. 

Finally, the theory is illustrated with an examination of the recent 
literature on lunar and interplanetary midcourse guidance. 

The first construction is carried out from the viewpoint of linear vector 
spaces, but only elementary principles from that subject are used (cf. ref. 6 ), 
A convenient notion in the second approach is the pseudo -inverse of a matrix 
(refs. 7 , 8 ). The two constructions give equivalent solutions for the guidance 
laws, although the application of the second construction will usually require 
less labor. 


SYMBOLS 


A(t) various matrices 

A(t 2 ,t!) transition matrix, relating the deviation state at t 2 to the 
deviation at ti 
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Aj(t 2 ,ti) 

^■z 

B(t) 

Bi(t 2 ,tx) , 
Ba(t 2 ,ti ) > B 


submatrices of A(t a ,tx), j = 1, . . . , 
azimuth angle 
various matrices 

matrices relating states in S(t 2 ) to states in S(tx) 


D 

Fi(t) 

G(t), Gi(t), 
G 2 (t) 

I 

N(t) 


n 

P 

RA 

R(t) 

r 

s(t) 

s 

t 

tjp 

u 

V(t), W(t) 
v 


declination angle 

matrix relating state deviations to deviations in orbital 
parameters 

guidance law matrix and its submatrices 
identity matrix 

state space; set of all deviations of position and velocity 
from a reference position and velocity at time, t 

unit vector normal to orbital plane 

m-tuple of independent orbital parameters (p 1 , . . . , p m ) 
right ascension 

set of all position deviations occupied by the satisfactory 
orbits at time, t 

position vector 

set of satisfactory deviation states 

element of S(t) 

time 

reference arrival time 
unit vector 

matrices appearing in construction of guidance laws 
velocity vector; also used as a unit vector 
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X(t) 

x(t) 

a 

• 

P 1 j P2 
8( ) 

M- 

9 

n 
c ) T 
c } 


vector appearing in construction of guidance law, otherwise 
used as a unit vector 

vehicle position and velocity, in Cartesian component, rela- 
tive to a central tody 

state vector; six Cartesian components of position and veloc- 
ity deviations from a reference orbit; sometimes written 



earth* s rotation rate 

specified values of parameters p 1 , . . . , p m of reference 
orbit 

constants relating entry flight time and entry range angle 
small deviation from reference value of ( ) 
gravitational constant of central body 
entry range angle 

matrices occurring in landing site guidance law 

vectors 

transpose 

set of objects defined by contents of { } 

DEFINITIONS 


It is useful to begin with definitions of the primitive notions underlying 
the theory. 

The satisfactory orbits are all those that meet the mission constraints on 
vehicle motion. The reference orbit is one member of the set of satisfactory 
orbits; the choice of reference orbit may be based on one or more factors such 
as fuel consumption, flight time, etc. 

Since the launch system attempts to inject the vehicle on the reference 
orbit, it is frequently adequate to assume that the orbits of interest in the 
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mid.ccu.rse guidance problem are sufficiently close to the reference orbit to 
allow the equations of motion to be linearized about the reference. With this 
assumption, the theory need consider only deviations from the reference state. 

The state deviation, x(t), refers to a column of the six Cartesian compo- 
nents of position and velocity deviations from the reference state at time, t. 
It will sometimes be replaced by the convenient form (br, Sv) , in which the 
position and velocity deviations are stated explicitly. There is no loss of 
generality in using position and velocity deviations to define the state since 
these are related to the deviations of any other six independent parameters 
by means of appropriate transformations. Further, this definition of the 
deviation state is especially suited to the development. 

The state deviations at two different times on an orbit are related by 

x(t 2 ) = A(t 2 ,ti)x(ti) (l) 

where A(t 2 ,ti) is the transition matrix obtained by linearizing the equations 
of motion about the reference orbit. Its columns are derivatives of the state 
coordinates at t 2 with respect to the state coordinates at ti evaluated on 
the reference orbit. 

Certain sets of states pertinent to the discussion are next defined. 

The set of all deviation states at time t forms a six -dimensional linear 
vector space W(t) called the state space. Each element of W(t), when added 
to the reference state at t, defines an orbit. The state space N(t) is 
everywhere equivalent to the set of all sextuples and, hence, includes some 
arbitrarily large deviations from the reference state. Although equation (l) 
inadequately represents the actual motion when the deviations are large, the 
theory is consistent based on the model of equation (l). In a practical con- 
text this means the linearized theory may be applied only in guidance problems 
in which equation (l) adequately describes the orbits that occur. Such a 
restriction is, of course, implicit in the use of linearized equations of 
motion. 

The set of satisfactory states S(t) at time t is the set of all those 
deviation states in W(t) which define satisfactory orbits. The guidance 
process is an attempt to correct an arbitrary state in W(t) to a state in 
S(t). The set of orbits defined by S(t) is fixed at all times; no new 
satisfactory orbits are added or others deleted as the mission proceeds. 

Lastly, the set R(t) is the set of position deviations occupied at time 
t by the satisfactory orbits. It is evident that any arbitrary state (Sr,Sv) 
can be corrected to a state in S(t) by means of an impulsive velocity cor- 
rection only if &r is in R(t) . 
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ANALYSIS 


Linear Impulsive Guidance Laws 

In midcourse guidance of space missions, it has frequently been assumed 
that the velocity corrections are impulsive, and that the orbits of interest 
can be described by linearized equations of motion. With these assumptions, 
the relationship between the state deviation and the velocity correction, 
that is, the guidance law, is a linear one and is defined over the entire 
state space. The class of guidance laws with these general properties may be 
termed linear impulsive guidance laws. 

One approach to the analysis is axiomatic; we begin with a definition of 
the class of linear impulsive guidance laws in terms of the elementary proper- 
ties it should possess, and then derive other properties of this class as a 
consequence of the definition. 


Definition . - A linear impulsive guidance law is a rule which specifies a 
unique instantaneous change of velocity, linearly related to the present state 
deviation, such that the final state is satisfactory. Such a correction is 
specified for every state; if the present state is satisfactory, then no cor- 
rection is made. 


Suppose such a guidance law is defined at time t, and x(t) is any state 
deviation. The change of state specified by the guidance law is linearly 
related to x(t) for every x(t): 


Ax = G(t)x(t) 

The 6x6 matrix G(t) will be termed the guidance law matrix at t. Since 
only the velocity is to be changed, then G(t) has the form 


G(t) 


0 0 
Gi(t) G 2 (t) 


where G(t) has been partitioned into 3X3 submatrices. The velocity correc- 
tion can also be written as 


Av = Gi(t)5r(t) + Ga(t)5v(t) 


After a correction has been made, the final state is 

x^(t) = [I + G(t)]x(t) 

If the present state is satisfactory, then no correction is required; that is, 

G(t)x(t) = 0 
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when x(t) is a satisfactory state. Since [I + G(t)]x(t) is a satisfactory 
state for any x(t), the guidance law matrix satisfies the equation 

G-(t) [I + G(t)] = 0 


The partitioned form of 
to 


G(t) can "be inserted to reduce the previous equation 


G-i(t) = -G2(t)Gi(t) 

^(t) = -G2 2 (t) 


( 2 ) 


The guidance law matrix associated with linear impulsive corrections satisfies 
equations (2) , hut not all matrices which satisfy equations (2) are necessarily 
guidance law matrices. The next section will determine existence conditions 
and the guidance law matrix in terms of various features of the set of satis- 
factory states. 


Existence Conditions and the General Guidance Law Matrix 

This section explores the connections between the guidance law matrix and 
the set of satisfactory states. It is found that linear impulsive guidance laws 
occur when the satisfactory states are certain types of subspaces of the state 
space , and that the guidance law matrix can be given in terms of appropriate 
states from S(t) . Some general properties of the guidance law matrix are 
given and, in particular, it is found that the submatrix Ga(t) can have only 
certain mathematical forms. 

The satisfactory states .- Suppose the guidance law exists and G(t) is 
its matrix. Then from the definition of the guidance law, the states in 
[I + G(t)]N(t) are satisfactory and also 

G(t)x(t) = 0 

for any state which is satisfactory. It follows easily that the converse of 
these two statements is also true. Consequently, if G(t) exists, then the 
satisfactory states can be defined from the guidance law matrix by either of 
the following 


S(t) = [I + G(t)]N(t) 

S(t) = (x(t) : G(t)x(t) = 0} 


( 3 ) 


Necessary con ditions for the existence of G(t) . - The first of equa- 
tion ( 3 ) states that if G( t ) exists, then S(t ) is the linear vector sub- 
space of N(t) formed by the set of all linear combinations of the columns of 
I + G(t ) ; that is, the columns of I + G(t) are satisfactory states and span 
the sub space S(t) . The dimension of S(t) is equal to the order of the 
largest square nonsingular submatrix of I + G(t) . In terms of its 3X3 
submatrices 


7 



I + G(t) = 


I 


0 


lGi(t) I + G 2 (t)J 

and it is clear that the dimension of S(t) is at least 3- 


(4) 


Further, R(t), the set of positions occupied at t by the satisfactory 
orbits is, equivalently, the set of all linear combinations of the position 
elements of any set of states which span S(t) . In particular, R(t) is given 
by all linear combinations of the columns in the upper half of I + G-(t) and 
is therefore three-dimensional, since the upper left submatrix I in equa- 
tion (4) is nonsingular. These statements are collected in the following 
necessary conditions for the existence of a linear impulsive guidance law 
matrix . 


Lemma 1 : If G(t) exists, then S(t) is a linear vector subspace 

of W(t) such that R(t) is three-dimensional. 

Thus, S(t) is a three-, four-, or five -dimensional subspace of N(t). 

The trivial case in which S(t) is six -dimensional need not be considered, 
since in this case every state in N(t) is satisfactory, and no correction is 
ever made. That R(t) must be three-dimensional is physically obvious, since 
it guarantees that at least one satisfactory state exists corresponding to 
every position deviation; and in any other case there will be position devia- 
tions for which no satisfactory state exists and for which a correction is 
therefore impossible. 


The guidance law matrix.- The necessary conditions for existence given by 
lemma 1 are also sufficient. The proof of this fundamental fact will follow 
from lemma 2 and theorem 1 below. Lemma 2 is a formality; theorem 1 shows that 
the guidance law matrix can be given from any basis of the satisfactory states 
and will be useful in the later text. 

Lemma 2 : Let S(t) be an n-dimensional subspace of N(t) and let 

Si = (ri,Vi), . . ., s n = (r n ,v n ) be any basis of S(t). Then 
R(t) is three-dimensional if and only if 3 < n < 6 and three of 
the vectors ri, • . . , r n are independent. 

A basis of a linear vector space is any set of linearly independent ele- 
ments in the space which will generate every element in the space . The number 
of elements in a basis is equal to the dimension of the space. Further, the 
space can be generated as the set of all linear combinations of the basis ele- 
ments. Thus, S(t) is the set of all linear combinations of si, . . ., s n , 
and R(t), the set of all positions occupied by the satisfactory states, is 
the set of all linear combinations of ri, . . ., r n . Hence, R(t) is three- 
dimensional if and only if n > 3 and three of the vectors ri, • . . , r n 
are independent. 



I 


Theorem 1 : If S(t) is an n-dimensional subspace of N(t), 

n = 3, 4 , and. so. = (rijYi), . . . , s n = (r n ,vn) form a basis of 
S(t) such that ri, r 2 , F3 are independent (renumber the states 
if necessary), then a linear impulsive guidance law exists and its 
matrix is given as follows: Let W, V be the 3 X 3 matrices whose 

columns are Fi, r2, ¥3, and Vi, V2, V3, respectively. Then 


where ( i ) 

(ii) 


Gi(t) = - 

G 2 (t) W _1 




-I 

for 

n = 3 

G 2 (t) = 1 

-I + 11U 

for 

n = h 


— — T _ — T 

-X + U1U1 + "U -2^-2 

for 

LT\ 

II 

3 


G(t) specifies the minimum velocity correction 

u is a unit vector in the direction of the (nonzero) 
vector .. 

v 4 - W r 4 


(iii) Ui, U2 are any orthonormal pair of vectors in the plane 
of the two (nonzero, independent) vectors 


Vj - VW -1 r j , j = 4, 5 


The proof of theorem 1 is carried out in appendix A. Note that the matrix 
¥ is nonsingular since its columns are assumed independent. The form of G(t) 
is not the most general linear impulsive guidance law matrix when S(t) is 
four- or five -dimensional but corresponds to the condition that the magnitude 
of the velocity correction be minimized. The general result is noted in 
appendix A but is not of interest here except in discussing the uniqueness of 
the guidance law matrix. 


In view of lemmas 1 and 2 , if G(t) exists, then the conditions of 
theorem 1 can always be satisfied; therefore, the guidance law matrix can 
always be given in one of the forms of that theorem. 


It should be noted in theorem 1 that guidance law matrices exhibit certain 
properties of form which depend only on the dimensions of S(t) . The submatrix 
G's(t) can have only certain mathematical forms. If we write 

wj = Vj - VW -1 rj, 3 < j < n 

then G2 ( t ) can be described as the reflection-projection matrix with pro- 
jection to the space orthogonal to the n -3 independent vectors Wj . An 
alternate view is that the mission constraints are independent of the wj 
components of the velocity deviation. Thus, G 2 (t)Sv cancels out the part of 
Sv which is orthogonal to the Wj and leaves unchanged the components of &v 
along Wj 
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Such a restriction on the form of G 2 (t) is expected in view of 
equations (2) . However, a larger class of matrices than the reflection- 
projection matrices given in theorem 1 will satisfy equations (2); the remaining 
matrices do not reflect the condition that the velocity correction is to he 
minimized. That is, since the mission constraints are independent of the wj 
components of the final velocity deviation, then these components can he 
changed arbitrarily without affecting the mission constraints, and the minimum 
correction is the one that makes no change in these components . 

Further, the suhmatrix Gi(t) is restricted to the extent that the range 
space of GxCt) is included in the range space of G 2 (t); that is, Gi(t)&r 
will lie in the space orthogonal to the wj . These properties are partly sum- 
marized hy the following singularities, which follow easily from theorem 1: 

T T 

u Gi(t) = u G 2 (t) = G 2 (t)u = Av . u = 0 for n = k 

_T _T _ 

UjGi(t) = UjG 2 (t) = G 2 (t)uj = Av . vj = 0 j = 1, 2 for n = 5 


Necessary and suff icient conditions for the existence of G(t).- Finally, 
lemma 2 and theorem 1 together state the converse of lemma 1; that is, if S(t) 
is a linear vector subspace of N(t) such that R(t) is three-dimensional, 
then G(t) exists. These results combine to give a necessary and sufficient 
condition for the existence of the guidance law. 

Theorem 2 : A linear- impulsive guidance law exists at t if and only 

if S(t) is a linear vector subspace of N(t) such that R(t) is 
three -dimensional . 

Uniqueness of th e guidance law matrix.- The uniqueness of the guidance law 
matrix corresponding to a given set of satisfactory states is investigated in 
appendix B. If G(t) is assumed to exist, the results state that G(t) is 
unique, provided S(t) is three-dimensional and is unique when S(t) is four- 
or five -dimensional if the magnitude of the velocity correction is minimized. 

To summarize the results of this section, a linear impulsive guidance law 
occurs when S(t) is any three-, four-, five-, or, trivially, six -dimensional 
sub space of N(t) such that there is at least one element of S(t) corre- 
sponding to every position deviation. When this condition is satisfied and the 
velocity correction is assumed to be minimized, then the guidance law matrix 
is given uniquely by theorem 1 from any basis of S(t) . 


Invariant Properties and Propagation of the Guidance Law Matrix 

Since the guidance law is of continuing interest throughout the mission, 
this section will consider the propagation of the guidance matrix. We first 
determine whether the properties discussed in the previous section are invariant 
with time. It is found that S(t) remains a linear vector space of invariant 
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dimension throughout the mission. Hence, the form of G 2 (t) is also invariant 
for all times when G(t) exists. Secondly, if G(t 0 ), the guidance law at the 
start of the mission, is assumed to exist, the relation between G(t) and 
G(t 0 ) can be determined for those times at which G(t) exists. 

Invariant properties . - The set of satisfactory orbits is fixed throughout 
the mission. If x(t Q ) is on a satisfactory orbit, then from equation (l) the 
corresponding satisfactory state on the same orbit at time t is 

x(t) = A(t,t 0 )x(t 0 ) 

If G(t 0 ) exists, then S(t Q ) is a linear vector space. By tailing the union 
over all states in S(t D ), we obtain 

S(t) = A(t,t 0 )S(t 0 ) ( 5 ) 

The state transition matrix A(t,t 0 ) is never singular (ref. l); 1 hence, 
if G(t Q ) exits, then from equation (5), S(t) is everywhere a linear vector 
space with the same dimension. Therefore, the form of G 2 (t) and other 
dimension-dependent properties of form noted in theorem 1 are also invariant 
wherever the guidance law exists. Further, since S(t) is everywhere a linear 
vector space then the existence conditions of theorem 2 are met if R(t) is 
three-dimensional. It should be noted that there may be one or more times 
during a mission at which this condition is not met. In summary, 

Theorem 3 « If G(t 0 ) exists, then (i) S(t) is everywhere a linear 
vector space with time -invariant dimension, and (ii) G(t) exists 
if and only if R(t) is a three-dimensional linear vector space. 

Existence criteria .- In the following, G(t 0 ) is assumed to exist, and 
the relation between G(t) and G(t 0 ) is obtained. It is first necessary to 
give the existence condition in terms of G(t Q ) and the transition matrix. 

The transition matrix may be partitioned into its 3 x 3 submatrices: 


A(t 2 ,ti) - 


"Ai(t 2 ,t 1 ) A 2 (t 2 ,ti)~‘) 

_A3(t 2 ,ti) A 4 (t 2 ,t!) 


( 6 ) 


The satisfactory states at t, from equations (3) and (5), can now be expressed 
as 

Bi(t,t 0 ) A 2 (t,t 0 )[l + G 2 (t D ) ] 

B 2 (t,to) A 4 (t,to)[l + G 2 (t Q )] 


S(t) = 


N(t 0 ) 


(7) 


^Ef A(t,t 0 ) were singular, distinct states would exist at t G which 
would generate identical states at t under free motion and, hence, distinct 
orbits which become identical orbits. 
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(8) 


The following notation has "been adopted in equation ( 7 ) 

Bi(t,t 0 ) = Ai(t,t 0 ) + A s (t,t 0 )Gi(t 0 ) 

B 2 (t,t 0 ) = Aa(t,t 0 ) + A 4 (t,t 0 )Gi(t 0 ) 

The positions occupied by the satisfactory orbits at t form the linear 
vector space given by 

R(t) = { Sr ( t ) = Bi(t,t 0 )Sr(t 0 ) + A 2 (t,t 0 )[l + G 2 (t 0 ) ]Sv(t G ) ; 

Sr(t 0 ) ,5v(t 0 ) arbitrary} ( 9 ) 

A necessary and sufficient condition for the existence of G(t) is that 
r( t ) be three-dimensional; that is, at least one satisfactory state exists 
corresponding to every position deviation at t. The matrix G 2 (t 0 ) can have 
only certain forms as given in theorem 1 of the preceding section. These forms 
depend only on the dimension of S(t) 



-I 


for 

n = 

3 

G 2 (t 0 ) - 1 

-I 

— T 

+ uu 

for 

n = 

4 


-I 

T T 

+ UxU-x + U 2 u 2 

for 

n = 

5 


where n is the dimension of S(t) and u, u 1; u 2 are the unit vectors of 
theorem 1. 

Existence criteria are obtained by substituting these forms of G 2 (t Q ) 
into equation ( 9 ) and determining the conditions for which R(t) is three- 
dimensional. Equation ( 9 ) gives R(t) as the set of all linear combinations 
of the six columns of Bx(t,t 0 ) and A 2 (t,t D )[l + G 2 (t 0 )]; hence, R(t) is 
three-dimensional when three of the six columns are independent. The results 
are readily obtained and are summarized by the following corollary. 

Corollary 1 ; If G(t 0 ) exists and dim[s(t)] = n, n = 3 , 4, 5 then 

G(t) exists if and only if one of the following is satisfied: 

(a) Bi(t,t 0 ) is nonsingular 

Cb) n = 4, Bi(t,t 0 ) has rank 2, and A 2 (t,t G )u is not in the 
range space of Bi(t,t 0 ) 

(c) n = Bi(t,t 0 ) has rank 2, and one of the vectors A 2 (t,t 0 )ui, 
A 2 (t,to)u 2 is not in the range space of Bx(t,t Q ) 

(a) n = 5 , Bi(t,t 0 ) has rank 1, and the two vectors, A 2 (t,t 0 )ux, 
A 2 (t,t 0 )u 2 are neither in the range space of Bi(t,t 0 ) nor 
dependent 



A propagation formula . - The guidance law matrix G-(t) exists as specified 
in corollary 1. When G(t) exists, it can he related to G(t 0 ) through the 
transition matrix. This relation is obtained below for case (a) of corollary 1 
in which B 1 (t,t 0 ) is nonsingular . The remaining cases are neglected here 
because the results have not been obtained in a useful form. 

Equation ( 7 ) gives S(t) as the set of all linear combinations of the 
columns of a matrix. To obtain the propagation formula, substitute the appro- 
priate form of G 2 (t 0 ) according to the dimension of S(t) into the matrix 
in equation ( 7 )- If n is the dimension of S(t), the resulting matrix must 
have n independent columns. For example, if n = 4 the matrix becomes 

— ij" 

Bi(t,t 0 ) A 2 (t,t D )uu 

_B 2 (t,t 0 ) A 4 (t,t 0 )uu T _ 

There is exactly one independent column among the right-hand three columns and 
this together with the left-hand three columns provides four independent states 
in S(t) . If, next, Bi(t,t 0 ) is assumed nonsingular then the conditions of 
theorem 1 are met; the matrix above provides four independent states in S(t) 
of which the first three provide three independent position deviations . Con- 
sequently, the guidance law matrix follows immediately from theorem 1. First, 
define the matrix 


B = A 4 (t,t Q ) - B 2 (t,t 0 )B 1 ;L (t,t 0 )A 2 (t,t 0 ) (lOa) 

This matrix is nonsingular (appendix C), and if A 2 (t,t Q ) is also nonsingular, 

B can be reduced to (appendix C) : 

B = [A 2 (t,t 0 ) _1 ] T B£ 1 (t,t 0 )A 2 (t,t 0 ) (10b) 

Finally, the guidance law matrix itself is given by the following corollary to 
theorems 1 and 3 '• 

Corollary 2 : If G(t Q ) exists and Bi(t,t c ) is nonsingular, then 

G(t) exists and is given as follows: 

Let S(t) be n -dimensional; n = 4, 5- Then 


Gi(t) - 

-G 2 (t)B 2 (t,t 0 )B£ 1 (t,t 0 ) 




’ -I 

for 

n = 3 

G 2 (t) = 

-I + u(t)iT I '(t) 

for 

n = 4 


-I + Ui(t)u! T (t) + u 2 (t)u 2 ’ I '(t) 

for 

n = 5 
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the form of G 2 (t) assumes that the velocity correction is 
minimized when n = k, 5 

For n = k, u(t) is a unit vector in the direction of 
Bu(t 0 ) 

For n = 5* Ui(t)j u 2 (t) are any orthonormal pair of 
vectors in the plane of the two (independent) 2 vectors 
Bux(t 0 ) , Bu 2 (t 0 ) 

u(t Q ), ui(t 0 ), u 2 (t 0 ) are the unit vectors associated 
with G 2 (t Q ) in theorem 1 and Bx(t,t 0 ), B 2 (t,t Q ), B 
are matrices defined in equations (8) and (lO) 


Derivation of Guidance Laws From Constraints 

The preceding analysis hegan with the primitive notion of satisfactory 
states and a definition of the linear impulsive guidance law and then obtained 
the theoretical results as consequences of the definition. In practice, guid- 
ance laws are obtained by applying constraints to the vehicle motion in order 
to meet mission requirements. This process is sufficiently general that the 
notion of guidance constraints must play an important role in the theory. 

In this section the connection of guidance constraints with the preceding 
theory is examined, and a uniform method for deriving linear impulsive guidance 
law matrices from the guidance constraints is given. 

Guidance constraints . - The guidance constraints are a number of restric- 
tions on the vehicle motion to be satisfied in order to meet the mission 
requirements. The guidance law is itself a statement of these constraints. 

In equation (3)> 

S(t) = (x(t) : G(t)x(t) = 0} (3) 

That is, those states which satisfy the mission requirements are given by equa- 
tion ( 3 ) as all states for which the linear function G(t)x(t) is zero. The 
number of independent scalar constraints imposed by the guidance law on the 
vehicle motion is then the number of linearly independent equations in 

G(t)x(t) = 0 

or, equivalently, the rank of G(t) . It is readily determined from equations ( 2 ) 
that the rank of G(t) is equal to the rank of G 2 (t), which, from theorem 1, 
is specified by the dimension of S(t) . 

Lemma. 3 : Let G(t) be the matrix of a linear impulsive guidance law 

defined at t. Then S(t) is an n-dimensional subspace of N(t), 
n = 3, k, 5, and the guidance law imposes 6-n independent scalar 
constraints on the state. 

2 Since Ux,u 2 are independent and B is nonsingular, Bux, Bu 2 are 
independent . 


where (i) 

(ii) 

(iii) 

(iv) 



It may "be noted that a linear impulsive guidance law cannot impose more 
than three constraints. Further, since n is invariant throughout the mission 
(theorem 3 ), the number of constraints imposed by G(t) is also invariant. 

Derivat ion of the guidance law from constraints ■ - In practice, guidance 
laws are obtained by applying constraints to the vehicle motion in order to 
satisfy the mission. This process is generally applicable and, with suitable 
restrictions, leads directly to linear impulsive guidance laws. 

Let 

p = (pi^ * • •} Pg) (n) 

be any six independent scalar parameters of orbital motion. This means that 
the set of all p, {p} , is in one-to-one correspondence with the set of all 
orbits . The six Cartesian coordinates of position and velocity at some time 
t, x(t), are often used, but, in general, the parameters may be chosen for 
convenience. In this case, X(t) can always be given in terms of p by virtue 
of their one-to-one correspondence with the set of all orbits: 

X(t) = F(p x , . . ., P s ; t) (12) 

The guidance constraints are a number of restrictions on the vehicle 
motion generated by the mission requirements. Suppose the guidance constraints 
have been expressed as constraints on independent orbital parameters and that 
these are included among the six parameters of (ll); in particular, suppose 
the mission requirements are met by all orbits for which the parameters 
p n+1 , . . ., p Q , n < 6, have the specified values a n+1 , • • cc.g. (if 
n =6, then no parameter is specified, and we have the trivial case in which 
all orbits satisfy the guidance constraints.) The satisfactory orbits are now 
defined by 

Ps = (V x , • • V n , cc n+1 , • • ., a s ); P x , • • •, P n arbitrary 

(13) 

A reference orbit is chosen from among the p s of equation (l3)> The remain- 
ing orbits can be referred to the reference orbit by: 

Sp = P - P ref = (8P 1# • • &P S ) (U) 

and the guidance constraints become 

&P n+1 = • • • = SP 6 = 0 (15) 

Assume, next, that the orbits of interest are sufficiently near the 
reference orbit that they may be adequately represented by the first two terms 
of a Taylor series expansion of equation (12) about the reference orbit. The 
deviation state at any time can then be expressed as: 
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x(t) = X(t) - X(t) 


ref 


Sr(t) 

SPi 

Sv( t) 

L 


Sr(t)~ 


Sp f 


e 

Sv(t) 

6p« 


Sp = F x (t)6p 


( 16 ) 


The matrix, Fi(t), is nonsingular by virtue of the one-to-one correspondence of 
N(t) and {Sp} with the set of all (linearized) orbits and, hence, with each 
other. That is, the columns of F x (t) are six independent states in N(t). 

The satisfactory states are all those deviation states that satisfy ( 15 ): 


S(t) =( 


x 


(t) = 


Sr(t) 

Sp 1 

Sv(t) 

Sp x 


Sr(t) 


6p n 

Sv(t) 


Sp 


n 



Sp 


i> ■ • 

arbitrary 


Bp. 


n 


(17) 


Equation ( 17 ) gives S(t) as the n-dimensional linear vector subspace of 
h(t) for which the n states 


Sr(t) 

Sv(t) 


Sr(t) 

Sv( t) 

L 

' a Pl J 

, . . . , s n - 

. 6p n 

^Pn _ 


( 18 ) 


are a basis. The existence conditions of theorem 2 for a linear impulsive 
guidance law can now be satisfied if and only if R(t) is three-dimensional. 
This occurs if and only if (i) n is not less than 3> and (ii) three of the 
vectors 


Sr (t) 


i- = 1 , 


n 


are independent. Assuming that these conditions are met, then the conditions 
of theorem 1 are fully satisfied by the n states of equation (l8) and the 
guidance law matrix can now be written directly from theorem 1. The following 
theorem gives the results : 

Theorem 4 : Let p^ = (p x , . . ., p s ) be six independent scalar 

parameters of orbital motion; let X 0 ('t) "be the state on a ref- 
erence orbit, and let the equations relating X(t) to p be 
linearized about the reference orbit: 


x(t) 


Sx(t) 

- 6p 1 


Sx(t) 

^6 


Bp 
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Then (i) a linear impulsive guidance law satisfying the constraints 

Sp j =0 j = n + 1, . . . , 6 

exists if and only if n > 3 and 


rank 


8r(t) 
- 5 Px 


dr(t) _ 

8p n J 


(ii) If the guidance law exists, then the guidance law matrix 
which specifies the minimum velocity correction is given 
as follows. Let Br/Bp^; k = 1, 2, 3 "be three independent 
vectors (renumber the parameters if necessary) and define 


then 


W(t) = 


Br( t ) Br(t) Br(t) 
5p 1 Bp 2 Bp 3 


V(t) = 


Bv(t) Bv(t) Bv(t) 


Bp, Bp, 


Bp, 


3 -I 


Gi(t) = -G 2 (t)V(t)w(t) _1 


G 2 ( t ) - ( 


-I for n 

_JP 

-I + uu for n 

T T 

-I + U]Ui + u 2 u 2 for n 


3 

4 

5 


where (i) all derivatives are evaluated on the reference orbit 

(ii) for n = 4, u is a unit vector in the direction of 
(Bv/Bp 4 ) - VW~ 1 (8r/8p 4 ) 

(iii) for n = 5* Ui> u 2 are any orthonormal pair of vectors 
in the plane of (Bv/Bp k ) - VW _1 (Br/Bp k ) ; k = 4, 5 


The results in theorem 4 parallel those of theorem 1 except that the guidance 
law matrix is given in terms of appropriate derivatives of the state. The 
derivatives appearing in equation (l8) and in the guidance law matrix form a 
basis of S(t), but the guidance law is unique and independent of the basis of 
S(t) that is used; that is, the derivatives appearing in the guidance law 
matrix can be taken with respect to any n independent scalar parameters of 
orbital motion which are also independent of the constraint parameters. 
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It may be noted that a linear impulsive guidance law does not occur if 
more than three constraints are required as, for example, in the rendezvous 
constraints [Sr(tf) = 0, Sv(tf) = 0], 


An Alternate Construction of Guidance Law Prom Constraints 

Thus far, linear impulsive guidance laws have been treated by considering 
the set of satisfactory states; theorem 1 shows that G(t) can be given from 
any basis of S(t) . To relate the guidance law to the guidance constraints, 
according to the previous section, if there are m (m < 3) independent guidance 
constraints, then a basis of S(t) is given by derivatives of the state with 
respect to any 6-m independent parameters which are also independent of the 
constraint parameters. The guidance law matrix can then be formulated from 
these derivatives . 

An alternate construction of the guidance law matrix from the guidance 
constraints is obtained by considering directly the conditions under which the 
constraints can be satisfied. The resulting definition of the guidance law 
will usually be simpler to apply than theorem 4. 

Let pT = (p , . . ., Pq) be m independent scalar parameters of orbital 
motion related to the present state X(t) by 

p = F[X(t) ] (19a) 

and suppose that the guidance constraints are that p is the specified vector 

P T = (am, . . • , %) 

Choose a reference orbit X Q (t) which satisfies these constraints and assume 
that the orbits of interest are sufficiently close to the reference orbit that 
the parameters p may be approximated by the first two terms of a Taylor 
series expansion of equation (l9a) about the reference orbit. Subtracting the 
reference equations from (l9a), obtain 

Sp = A(t)Sr(t) + B(t)8v(t) (l9b) 

Here, A(t), B(t) are m X 3 matrices in which the jth rows are the deriva- 
tives, evaluated at the reference orbit, of pj with respect to the Cartesian 
coordinates of position and velocity, respectively. Since the m parameters 
of p are independent, then 

rank [A B] = m (20) 

The guidance constraints are 

&Pj = 0 j = 1, • • -, m (21) 
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and we define the existence of a guidance correction which satisfies these 
constraints as follows : 

A guidance correction, Av, exists at t if for every deviation 
state (5r, Sv) there exists a Av which satisfies the m 
independent equations : 

0 = A(t)5r + B ( t ) ( 5 v + Av) (22) 

Exist ence conditions . - An existence condition for the guidance correction 
is readily established. 

Theorem 5 : Equation (22) has a solution (i.e., a guidance correction 

exists at t) if and only if 

rank (B) = m < 3 

Proof - Equation (22) can be rewritten as: 

B(t) Av = -[A(t) B(t)]Sx(t) (23) 

where A, B are m X 3 matrices. We seek the conditions for which these m 
inhomogeneous linear equations have a solution, Av, for arbitrary 6x. The 
Krone cker-Capelli theorem (ref. 6) states that equation ( 23 ) is compatible 
(has at least one solution) if and only if the rank of B equals the rank of 
the augmented matrix, B a , where the augmented matrix is formed by attaching 
the m x 1 column 

b = [A B]6 x 

to the coefficient matrix, B; 

B a = [B b] 

Thus, the Kronecker-Capelli theorem is satisfied if and only if b is a linear 
combination of the columns of B. But b can be any linear combination of the 
columns of [A B] since 8x is arbitrary; hence, the columns of [A B] are 
required to be linear combinations of the columns of B. Thus, 

rank (b) > rank [A B] = m 

and, since the columns of B are contained in the columns of [A B], 

rank (B) < rank [A B] 

Thus, the Kronecker-Capelli theorem is satisfied if and only if 

rank (B) = rank [A B] = m (24) 

Note that B is an m X 3 matrix whose maximum rank is minimum fm,3} and 
equation (24) cannot be satisfied if m > 3* 

T he guidance law matrix .- Suppose that the guidance correction exists. 

Then the general solution of equation (22) can be given directly from refer- 
ence 7 (theorem 2, corollary l) ; 

Av = -B + ASr - B + BSv + z (25a) 
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where B+ indicates the pseudo -inverse 3 of B and the vector z is any 
arbitrary vector orthogonal to the domain (row space) of B. Now, the first 
two terms on the right-hand side of equation (25a) lie in the domain of B; 
in fact, B+B is the matrix which projects on to the (m-dimensional) domain of 
B. Thus, z is orthogonal to the first two terms, and if we adopt the condi- 
tion that z is to he chosen such that Av is a minimum, then z will he 
zero. The desired solution of equation ( 23 ) is then 

Av - -B + A6r - B + B5v ( 25 b) 

The corresponding guidance law matrix can now he identified: 

Gi(t) = -B + A 

0 2 (t) = -B + B 

It is readily verified that equation ( 25 b) satisfies the definition of a linear 
impulsive guidance law. Further, that definition can he satisfied by more 
general choices of z than z = 0, as noted in appendix A, hut the general case 
is not of interest in this paper. 

Thus far, it has been shown that if a guidance correction exists , a matrix 
can he written for a linear impulsive guidance law. Conversely, if that matrix 
can he written, then equation ( 23 ) has at least one solution for every devia- 
tion state, Sx, and the guidance correction can he said to exist. In summary, 

Theorem 6 : Let pT = (p 1 , . . ., p^ he m independent scalar 

parameters of orbital motion; let X Q (t) he the state on a reference 
orbit, and let the equations relating p to the state he linearized 
about the reference orbit 

6p = A(t)Sr(t) + B(t)Sv(t) 

Then (i) a linear impulsive guidance law satisfying the constraints 
Sp = 0 exists if and only if 

rank B(t) = m < 3 

(ii) If a guidance law exists, then the guidance law matrix which 
specifies the minimum velocity correction is 


3 The pseudo -inverse is a generalization to arbitrary matrices of the 
inverse of nonsingular square matrices. References 7 and 8 may he consulted 
and some pertinent examples of pseudo-inverses are given in the next section 
of this paper. The pseudo-inverse, B + , of B is the unique solution of the 
four equations : 

BB+B = B, BB + = (BB+) 1 
B+BB+ = B+, B+B = (B+B) T 

Note that if B is a nonsingular square matrix, then B+ is the ordinary 
inverse, B _1 . 
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Gi(t) = -B + (t)A(t) 

G 2 (t) = -B + (t)B(t) 

Theorem 6 will usually be easier to apply than theorem 4. It should be 
noted that only the 6m (scalar) derivatives of the constraint parameters with 
respect to the Cartesian coordinates of the state are required for the con- 
struction of G(t) • Theorem 4 requires the completion of the m constraint 
parameters to a set of six independent parameters and then defines G(t) from 
6(6 -m) derivatives of the state coordinates with respect to the new 
parameters. 

It may be noted that G(t) can be expressed in terms of the matrices 
A(t 0 ), B(t 0 ) for some fixed time, t Q , and the transition matrix, A(t,t 0 ) from 
the following relations 

A(t) = A(t 0 )A!J - B(t 0 )Al 
B(t) = B(t 0 )A? - A(t 0 )A 2 

where Ax, . . . , A 4 are the submatrices of A(t,t 0 ) as in equation (6). 

Remarks_ on theorems 4 and 6 . - Theorems 4 and 6 are necessarily equivalent; 
it can be demonstrated independently that the existence conditions and guidance 
law matrices of these two theorems are equivalent. The method is outlined 
briefly below. 


Suppose p x , . . . , p s are six independent scalar parameters of orbital 
motion related to the state deviation by the linearized equations : 


5p = A(t)Sx(t) 
6x(t) = B ( t ) Sp 


(26) 


where A, B are nonsingular 6x6 matrices and B = A -1 . Suppose that the 
guidance constraints are 


and partition A, B 



11 

•i"D 

ft 

cO 

0 

j = 

= 1, • • 

m < 3 


as follows : 






Ax 

a 2 


(m X 3) 

(m X 

3) 

A = 

_A3 

a 4 


_(6-m X 3) 

(6-m 

x 3)_ 


— 

— 


_ 


— 


Bx 

b 2 


(3 x m) 

(6-m 

x 3) 

B = 

_b 3 

b 4 _ 

— 

_(3 X m) 

(6-m 

x 3). 


Theorem 6 constructs the guidance law from Ax and A 2 , and theorem 4 gives 
G(t) from B 2 and B 4 . Noting that AB = BA = I e , we have the two relations 


21 



AiB 2 + A 2 B 4 = 0 (m x 6-m) 

BiA 2 + b 2 A4 = 0 (3x3) 


(27a) 

(27b) 


With these two relations it can "be shown that (proof omitted) for m < 3: 

rank (A 2 ) = m IFF rank (B 2 ) = 3 

The left- and right-hand sides above are, respectively, the existence condi- 
tions of theorems 6 and 4; hence, the two are equivalent. 

Next, assume that the existence conditions are satisfied and let W he 
the matrix whose columns are the three independent columns of B 2 and V the 
corresponding columns of B4. Further, let r 1 , Vj he any remaining column 
of B 2 and its corresponding column from B 4 . Then equation (27a) yields 

-A 2 + Ax = (a 2 + a 2 )w _1 

-Ae + Airj = (A 2 + A 2 )vj 

Substituting (28a) into the second equation yields 

[A 2 + A 2 ] [ vj - W -1 rj ] = 

From theorem 6, G 2 (t) is the m-dimensional reflection-projection matrix pro- 
jecting to the space orthogonal to the 3~m independent vectors, Vj - W -1 rj. 
From equation (28h) it readily follows that the result for G 2 (t) is identical 
in the two theorems. From this result and equations (28a) it follows immedi- 
ately that the expressions for Gx(t) are also identical in the two theorems. 


(28a) 

(28b) 


Examples 

Some examples of guidance law matrices from the lunar and interplanetary 
mission studies are given below. These will serve to illustrate the applica- 
tions of theorems 4 and 6 and the properties of form for matrices of linear 
impulsive guidance laws. 

Example 1; Fixed time of arrival guidance (e . g. , ref, l) . - The guidance 
constraints require an interception of the reference orbit at some fixed t ime , 
tp. The linearized relation between the state at two different t ime s is given 
by equation (l). With the transition matrix partitioned as in equation (6), the 
position deviation at tp is: 

5r(tp) = Ax(tp,t)Sr(t) + A 2 (tp,t)Sv(t) (29) 

The linearized guidance constraints are: 

Sp = 8r(tp) = 0 
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Theorem 6 can be applied directly to obtain: 

(a) G(t) exists if and only if As(tp,t) is nonsingular 4 

(b) If G-(t) exists, then 

Gi(t) = -Ai 1 (t F ,t)A 1 (t F ,t) 

G 2 (t) = -I 

This is the principal example in the class for which S(h) is three- 
dimensional. Other examples can be given by applying any three independent 
scalar constraints to the vehicle motion. 

A number of variable arrival time guidance laws occur to illustrate the 
class for which S(t) is four -dimensional. Some of these are based on the 
principle of nulling only part of the miss (position deviation at arrival) 
rather than using fixed-time -of -arrival guidance to correct the total miss 
(e.g., refs. 2 , 2 >, 4). If such a guidance law can be used, the fuel require- 
ments of fixed-time -of -arrival guidance will usually be relieved, as in some 
applications of reference 2 which have employed a single midcourse correction, 
or in the simulations of reference 3 where the downrange miss cannot be esti- 
mated very well until late in the mission. 

Example 2 : Interplanetary variable arrival time guidance (ref. 4) . - The 

guidance attempts to obtain a desired position relative_to the target planet. 
The reference orbit has the desired relative position, A c , at the reference 
arrival time tp. Let V R be the velocity of the vehicle relative to the 
planet at tp. Then the relative position of the vehicle at times near tp 
on any nearby orbit can be approximated by the (linearized) expression 

S(tp + St) = Ao + Sr(tp) + V R St (30) 

This expression assumes that the relative velocity is approximately fixed on 
all nearby orbits (acceleration of the vehicle and planet in the time interval 
(%> ,t F + St) has been neglected) . The constraints are satisfied if Br(tp) 
is chosen so that 

A(tp + 5t A ) = Aq 

for some delay in the arrival time, bt A . For all orbits on which this is true, 
equation ( 30 ) gives the relation 


Sr(t F ) = -V R &t A ( 31 a) 

that is, the constraints are satisfied by all orbits for which Br(tp) is 
parallel to V R . The converse is also true. Equivalently, the guidance seeks 

“^The matrix As(tp,t) is singular in Keplerian orbital motion when t 
and tp ' either are identical or correspond to points of the reference orbit 
which differ in true anomaly by rat. Such instances can occur in interplane- 
tary flights, as in the return legs of the orbits of reference 3« 
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to cancel the lateral miss (position deviation at t F perpendicular to Vr), 
hut the downrange miss (along Vr) can he arbitrary since it only determines 
the delay in arrival time. 

Let v he a unit vector in the direction of Vr; then the guidance con- 
straints can he written as 

Sp = (i - w^)Sr(tjp) = 0 (31b) 

Equation (29) gives the relation between the guidance constraints and the 
state deviations : 

(I - vv r )Sr(t F ) = (I - vv [ ')[Ai(t F ^t)&r(t) + Ae(tp,t)&v(t) ] 

Theorem 6 can he applied directly to give: 

(a) A guidance law exists IFF rank (i - vvT)A 2 (tp,t) = 2 
(h) If a guidance law exists, then (time arguments dropped): 

Gi(t) = -[(i - W r )A 2 ] + (l - W r )A i 
G 2 (t) = -[(i - w T )A 2 ] + (l - w T )A 2 

If A 2 is nonsingular, then (a) is automatically satisfied and the expressions 
in (h) can he reduced as follows: 

[(I - W r )A 2 ] + = (i - ww r )A 2 1 

where 

A 2 v 

w = 

jA^vl 

(i - wwl)A 2 1 (l - vv®) = (i - w^)A 2 1 


from which the guidance law matrix becomes 

Gi(t) = -(i - ww^ 1 ) A 2 1 (tj>,t )Ai(tp,t ) 
G 2 (t) = -I + ww 1 


(32) 


This result is similar to the fixed-arrival -time guidance except that the w 
component of velocity is unchanged since it affects only terminal position 
deviations along Vr. Equations (32) preserve the downrange miss due to 
injection errors, and if this miss is large, there will he large deviations 
from the reference orbit near the end of the mission. 



Example 3 : Term inal ph ase pericenter guidance (ref. 3 ) 5 - - Pericenter 

guidance Vas derived for use in the terminal phase of interplanetary missions. 
It operates within the sphere of influence of the target planet and seeks to 
obtain the reference value of the radius of pericenter and to place the peri- 
center position in the plane of the reference orbit relative to the targe + 
planet. Reference 3 derives the guidance law matrix from derivatives of the 
two-body equations of motion. An alternate formulation is given here in terms 
of the transition matrix. 

The linearized equations governing the pericenter position are of the form 

Sr p = R(t)6r(t) + V(t)Sv(t) (33) 

where (&r,8v) is the state deviation from the reference orbit relative to the 
target planet. In particular, if tp is the reference arrival time corre- 
sponding to pericenter on the reference orbit, then write equation (33) as: 

8r p = R 0 Br(tp) + V 0 8v(tp) ( 3*0 

where R 0 , V 0 relate Brp to deviations from the reference orbit at tp. 

The matrices R 0 , V Q are readily determined as follows. Define the 
orthonormal vectors u, v, n which lie, respectively, along the reference 
pericenter position, along the velocity vector at pericenter, and along the 
normal to reference orbital plane. The time t F corresponds to pericenter on 
the reference orbit so that the reference state at t-p can be written 

r(t F ) = r pQ U 

V(tp) = V Po v 

Next, consider, successively, small deviations Br(tp), Bv(tp) along the three 
directions u, v, n and obtain the corresponding deviations in pericenter 
position, Brp. Summing these deviations (algebra omitted), obtain 

Brp = (uu T + nn T )6r(t F ) - ( ^ ^vu 1 ^) Bv(tp) ( 35 ) 

The guidance constraints require that the pericenter position lie in the 
reference orbital plane, or 

-T r — n 
n Brp = 0 

and that the value of the radius of pericenter be the same as the reference 
value ; 

|r(t F ) + Br p | - |r(t F )| = 0 


5 The term, per lapse guidance, used in reference 3 , is lexicographically 
doubtful and has been changed to pericenter guidance in this paper. 
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This can he linearized to the form: 

_T _ 

u x Brp = 0 

Inserting these in equation (35 ) > the linearized guidance constraints are: 

-T 1 


U ~ ' Sr p = ' ^ 


Sr(t F ) = 0 


(36) 


-T , 
n ‘ 


n^ 1 


Equation (36) states that the constraints for pericenter guidance are 
equivalent to cancelling out only the lateral miss at t F . These are exactly 
the constraints for example 2 , and the resulting guidance law should therefore 
have the same form as equation (32) . 

Equation (29) may he substituted for Sr(t F ) in equation (36), and the 
linearized equations for the constraint parameters become : 



[Ai(t F ,t)Br(t) + A 2 (t F ,t)Bv(t) ] 


Theorem 6 can now be applied to obtain the guidance law matrix: 


— T ' 
u 


— T 

n 


A 2 ( bp > t ) 


= 2 


(a) A guidance law exists IFF rank 

(b) If a guidance law exists , then (time arguments dropped) , 

G-i(t) = - I I I Ao II I An 


— . 1 - 

+ / m \ 

/ ■ — T 1 \ 

/ -T\ 

/ U 1 

/ u 1 

1 1 -^-2 

1 

V — T/ 

\ — T / 

\ n / ^ 

\ n / 


G\a(t) - 


r_T\ 

u ' 


— T 

in 1 


-1 + 


Ac 


,'-T\ 
1 u 




Ac 


For this case A 2 is generally nonsingular, 6 whence (a) is generally satisfied 
and the expressions In (b) can be reduced as follows: 


■ 

fu 1 ! 



j?j 

a 2 

_ 

In / 

_ 


-1 + 


= (i - w^-)A 2 1 (un) 


s For Keplerian motion, the orbit relative to the target planet is hyper- 
bolic in which case A 2 (t F ,t) will be nonsingular except at t = tp. 
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where 


w = 


Aj T? 

i Ai^v 



From these, the guidance law matrix "becomes, as in equations (32); 

G i(t) = -(I - ww^A^Ai 
G 2(t) = -I + 


( 37 ) 


Exampl e 4 : Fixed- lan di ng- site guidance (ref. |?) «- This guidance law has 

"been applied in studies of the return phase of the lunar mission. The primary 
objective of guidance in this phase is to make a safe entry and, if possible, 
to land at a fixed landing site. Fixed-time -of -arrival guidance had been used 
in previous studies, with interception of the reference orbit at time tp 
corresponding to vacuum perigee on the reference orbit. Although adequate for 
the range of initial errors studied in reference 5, the FTA guidance law does 
not explicitly recognize the objectives of the midcourse guidance and for com- 
parison the fixed -landing -site guidance law was derived to satisfy the 
constraints : 


(a) Vacuum perigee altitude is fixed at the value of the center of the 
safe entry corridor. 

(b) The vehicle lands at a geographically fixed landing site without 
crossrange maneuvering during the entry flight. 


The second constraint is adopted in view of the limited crossrange maneuvering 
capability of Apollo-class vehicles, compared to the downrange maneuver 
capability. 


The guidance law matrix will be determined by defining the satisfactory 
deviation states at any time t and applying theorem 1. The state deviations 
at t are defined by their corresponding deviations at tp from equation (l): 

Sx(t) = A(t,tp)&x(tp) (l) 

The deviations at tp, in turn, can be defined by deviations in vacuum perigee 
position, velocity, and time; that is. 



To a first-order approximation 


6tp — tp - t^ 

r(tp) £ r(t p ) + (tp - tp)v(t p ) s F Pq + 5r p - Vp Q &t p 
V(tp) S V(tp) + (tp - t p )v(t p ) = Vp Q + 5V p + F p 5t p 

where the subscript, o, refers to values on the reference orbit. From these 
results, we have the 6 x 7 matrix: 


H = 


Is 





(39) 


The previous equations have transferred the problem of defining the satis- 
factory deviation states, Sx s (t), into one of finding perigee deviations, &r p , 
6v p , Btp, which satisfy the constraints. These are obtained as follows: Let 

u, v, n be unit vectors in the direction of vacuum perigee position, velocity, 
and the normal to the orbital plane. Then 


r p = r p u 


V p = V p v 


(40a) 


and u, v, n can be given in terms of the landing site location at the time of 
landing, and entry flight-path parameters by the following transformations to 
inertial coordinates: 


(uvn ) = (u L M) 


cos cp 


sin cp 


0 


-sin cp sin A Z p cos cp sin A 7iTi -cos A, 


■ZL 


-sin cp cos A y Ti cos cp cos A 7Ti -sin A Z p 


(u l M) = (IJK) 


cos Dp sin RAp -sin RAp 
cos Dp cos RAp cos RAp 
sin Dp 0 


-sin Dp sin RAp 
-sin Dp sin RAp 
cos Dp 


(40b) 


( 40c) 
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I 


The parameters are illustrated_in_ 
sketch ( a) . The vectors u L , E, N are 
the directions of the landing site and 
local East and North at the landing 
site; their components in the inertial 
frame are given "by equation (40c) and 
are functions of the landing site lati- 
tude, Dl, and right ascension, RA^Ct^) • 
For a geographically fixed landing site, 
the second transformation is a function 
of the time of landing only: 

RA L (t L ) = RA L (t Lo ) + a(t L - t Lo ) 

( 40d) 

where a is the earth's rotation rate. 


The remaining parameters, cp, Agp are the total range angle from vacuum 
perigee to landing, and the azimuth of the plane of the entry flight at the 
landing site latitude. Note that in the first transformation the entry flight 
motion is assumed to take place in a single inertial plane and, therefore, 
already contains the constraint that no crossrange maneuvering occurs during 
the entry flight. 

For all entry flights with entry speed near the escape speed, an approxi- 
mate linear relation between flight time and total range angle is available: 

t L “ tp = Pl<P + (^ 0e ) 

where (3 1 , P 2 are constants (ref. 5) • 

The satisfactory perigee conditions are found from equations (40a) to 
(40e) by imposing the constraints that rp is fixed at the value corresponding 
to the center of the entry corridor, rp Q , and that Dp, RApCtp) correspond to 
the desired landing site. The remaining parameters which may be varied arbi- 
trarily to define the satisfactory perigee conditions are then Vp , cp, Azl> 
and t]^. Assuming small deviations in these parameters about their values on 
the reference orbit, we obtain the satisfactory deviations in vacuum perigee 
conditions from the appropriate derivatives of equations (40a) to (40e) 
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where 


o' 

li> 

$ 

1 

r Po 

sin cp Q n 

ar-n Kxu 
-^0 


V 

V Po u 

- V Po 

cos cp Q n 

aV-n Kxv 
^0 

(41b) 

_o 

-Pi 

0 


1 



where u, v, n are now unit vectors taken on the reference orbit and the 
symbol, x, indicates the vector cross product. Finally, equations (l), (32), 
(41) are combined to give the satisfactory deviation states at t: 


5x s (t) = A(t,tpH 


I Bv p \ 

Sep 


SA, 


ZL 


iBtr 


(42) 


where 

* = 

Since A(t,tjp) is a nonsingular transformation of the state space, then the 
dimension of S(t) is the same as the rank of \|r, which can be shown to be 
four. Hence, the four-columns of A(t,tj>)t form a basis of S(t) and theo- 
rem 1 can be applied directly. Partition A(t,tp) as in equation ( 6 ); then 

(i) Cf(t) exists IFF rank [Ai A 2 ]\|r = 3 (43a) 

Assuming that (43a) is satisfied, let W be any three independent columns of 
[Ai A 2 ]f and V the corresponding columns of [A 3 A 4 ]i|/. Further, let ? 4 , V 4 
be the remaining column of these two matrices, respectively. Then 

(ii) G x (t) = -G 2 (t)W _1 

<J 2 (t) = -I 4- 5u T (43b) 

where u is a unit vector in the direction of V 4 - VW ^ 4 . 

These results are not in a satisfactory final form, but further discussion 
is beyond the scope of this paper. 

The guidance laws in the class for which S(t) is five -dimensional result 
from imposing a single scalar constraint on the vehicle motion. For instance, 
in example if the constraint is reduced to obtaining the desired pericenter 
radius without any adjustment of the plane of motion, the linearized guidance 
constraints become 


o = u^Srp = u^Ai(tp,t)&r(t) + u' I 'A 2 (tjn,t)Sv(t) 
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for which the guidance law matrix becomes (algebra omitted): 


Gi(t) - - 


— T 

su A3.(t F ,t) 
I Ap^ul 


--T 


Ga(t) = -ss 


where 


A 2 T (t ; g 1 ,t)u 

|A 2 '%| 


After the first correction, the corrected orbit tends to diverge from the 
reference orbit because the deviation of the actual orbital plane from the 
reference plane of motion is preserved; consequently, the linearized equation 
for u^5Fp becomes a poor approximation later in the flight. This difficulty 
is relieved when equations (37) define the guidance law matrix. 


CONCLUDING REMARKS 


This paper has investigated relationships among some of the basic notions 
of the linearized theory of midcourse guidance in space missions; the satis- 
factory orbits, linear impulsive guidance laws, and guidance constraints. 

A satisfactory orbit is one that satisfies the mission constraints on 
orbital motion; the set of such orbits is fixed for a mission, and the corres- 
ponding set of states at any time, t, is a subset of the state space called the 
satisfactory states. The guidance process is one of correcting the current 
state of the vehicle to a satisfactory state by means of a change of velocity. 

The guidance law is the rule that specifies the velocity correction, Ax, 
that should be made in order to correct the vehicle state, X(t), to a satis- 
factory state : 


Ax = G(x(t)) 

A simple class of such laws,- termed linear impulsive guidance laws, was studied 
in this paper. These guidance laws yield velocity corrections that are linear 
functions of the state deviation from the reference state, and are defined. on 
the entire state space : 


Ax = G(t)x(t) 

where G(t) is the matrix of the guidance law. In the study of midcourse 
guidance in space missions linear impulsive guidance laws arise when the mis- 
sion imposes three or fewer constraints on the vehicle orbital motion, and when 
the assumptions are made that velocity corrections are impulsive and the orbits 
of interest are sufficiently close to a reference satisfactory orbit to allow 
the equations of motion to be linearized about the reference orbit. 
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It was demonstrated that linear impulsive guidance laws can he defined 
from the satisfactory states (the converse is also true) . Linear impulsive 
guidance laws occur whenever S(t) is a three-, four-, five-, or, trivially, 
six-dimensional linear vector subspace of the deviation state space such that 
there is at least one satisfactory state corresponding to every position devia- 
tion (theorem 2) . In this case the guidance law matrix can he given from any 
set of basis elements of the satisfactory states (theorem l) . The guidance law 
matrix corresponding to a given set of satisfactory states is unique when S(t) 
is three-dimensional, and otherwise becomes unique if the minimum velocity cor- 
rection is specified. Theorem 1 shows that the guidance law matrix has several 
properties of form which depend only on the dimension of S(t); if S(t) is 
n- dimensional, then Ga(t) is a ref lection -pro jection matrix projecting to a 
(6-n) -dimensional space which is also the range space of Gx(t) . 


Some time -dependent properties of the guidance law matrix were considered. 
Since the set of satisfactory orbits is fixed, the satisfactory states at two 
different times in a mission are related by a nonsingular transformation 
(eq. ( 5 ))- Consequently, S(t) has a time invariant dimension, and for those 
times at which the guidance law exists, the dimension dependent properties of 
form of G(t) are also invariant. A propagation formula relating the guidance 
law matrix at two different times in the mission was also derived. 


It was demonstrated that the satisfactory states and the guidance law can 
be defined from the guidance constraints. In practice, a mission requires that 
the vehicle motion satisfy various constraints. If these are expressed in 
terms of independent orbital parameters 


Pj = ay j = 1, • • m 

where ctj are the specified values and m is the number of constraints, then 
the satisfactory orbits are defined as all orbits that satisfy these con- 
straints. If [pi, ■ . PjjJ are completed to a set of six independent 
orbital parameters and a satisfactory reference orbit is selected, then a basis 
of the satisfactory states, S(t), is given by the derivatives: 


Sx(t) 

"^ 7 ' 


j 


m + 1, 


6 


where the derivatives are evaluated on the reference orbit. Thus, theorem 1 
can be applied to define the existence conditions and the guidance law matrix 
in terms of this basis (theorem 4) . 


Finally, an alternate construction of the guidance law matrix is given 
from the constraints by considering the conditions for which a velocity correc- 
tion will satisfy the constraints. If the constraint parameters are expressed 
in terms of m independent orbital parameters and a reference orbit is 
selected, then for all orbits near the reference orbit, the deviations of the 
constraint parameters are related to the state deviations by the (linearized) 
equation: 

Sp = A(t)Sr(t) + B(t)&v(t) 


32 


II 



A correction is said to exist if for every state, Sx, a velocity change, Av, 
can "be found to satisfy the constraint equations: 

0 = A(t)Sr + B(t)(&r + Av) 

These are m-inhomogeneous linear equations for which an existence condition 
(theorem 5) and the general solution for Av (eq. (24)) are readily given. A 
linear impulsive guidance law matrix can then "be given from the general solu- 
tion (theorem 6) . 

The constructions of theorems 4 and 6 are different; theorem 4 gives G-(t) 
in terms of the derivatives, Sx(t)/Sp^, while theorem 6 uses the derivatives 
Sp-^/dx and theorem 6 will usually "be easier to apply in specific midcourse 
guidance problems. Nevertheless, the two results are equivalent. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Dec. 8, 1965 
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APPENDIX A 


PROOF OF THEOREM 1 


The results quoted in theorem 1 are included in the proof of the following 
abbreviated statement. 

Theorem : If S(t) is an n -dimensional sub space of N(t), 

n = 3 , 5 , 6, and Si = (ri,Vi), . . .,_ s n j= (r n ,v n ) are n 

independent states in S(t) such that ri, r 2 , r 3 are independent, 
then a linear impulsive guidance law exists. 

By hypothesis, S(t) can be a three-, four-, five-, or, trivially, six- 
dimensional subspace of N(t) and these cases will be treated separately 
below. 


Case a: n = 3 


Let Si, S2, s 3 be three independent states in S(t) 

si = (ri,vi), s 2 => (r 2 ,v 2 ), s 3 = (r 3 ,v 3 ) 

and assume ri, r 2 , F3 are independent vectors. Next, define W and V to be 
3X3 matrices whose columns are Fi, r 2 , F3 and Vi, v 2 , v 3 , respectively; note 
that W is nonsingular since its columns are independent. 


Let (Sr,Sv) be any state deviation. To correct this state to a satis- 
factory state first determine the satisfactory state (s) which has the position 
deviation Sr; that is, determine all those states (Sr,Sv s ) that are linear 
combinations of the form 



CCj.S 2 _ + + CO3S3 


(Al) 


where the scalars aq are chosen such that 


Sr = airi + a 2 r 2 + 0,3^ = W 



(A2) 


The constants of combination have a unique solution which is obtained by 
inverting equation (A2) : 

l ai \ 

a 2 = W _1 5r 

\a 3 / 

The corresponding velocity deviation on the satisfactory orbit which passes 
through Sr is, from equation (Al) : 


3 ^ 



or 


(A3) 


SVg = CXi"Vi + <X 2 "V 2 + a 3^3 

Bv s = VW -1 &r 

The vector, Sv s , depends uniquely on 8 r, and there is exactly one satisfactory 
state (&r,Bv s ) corresponding to each position deviation, 5r. The union of 
all such states is : 

S(t) = [x(t) = (&r,VW - 1 Sr); &r arbitrary] (A4) 

the velocity change required to correct the arbitrary state (Br,Sv) to the 
satisfactory state ( 8 r, 8 v s ) is then 

Av = Bv s - Bv = W 1 Br - Bv (A5) 

and the guidance law matrix can now be written as : 


Gi(t) = VW " 1 
G 2 (t) = -I 


(A 6 ) 


The guidance law obtained specifies a change of state to a satisfactory state 
for every arbitrary state in N(t) by means of a unique impulsive velocity 
correction linearly related to the original state. Further, no correction is 
made if the original state is satisfactory. In short, equation (A5) satisfies 
the definition of a linear impulsive guidance law. 


Case b : n = 4 

Let S]_, s 2 , s 3 , s 4 be any basis of S(t) and assume that ri, r 2 , r 3 
are independent. Define the 3x3 matrices, V, W as above. 

Next, let (&r,Bv) be any state in N(t) and determine all satisfactory 
states corresponding to 8r; that is, determine all linear combinations 


'Sr 


^5v c 




such that 


Sr = ocari + a 2 r 2 + a 3 r 3 + 0^4 


or 


rai\ 

Sr = W ( a 2 I + a 4 r 4 

\CL 3 J 


(AT) 


(AB) 
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The solution of (AS) is not unique, "but ax, a 2 , a 3 can he given in terms of 
the arbitrary parameter, a 4 , by inverting equation (A 8 ): 

I a 2 I = W -1 Sr - a A W - 1 r 4 ; a arbitrary 

Vs/ 

The corresponding velocity deviations on the satisfactory orbits which pass 
through Sr are, from equation (A7): 

Sv s = cxxVx + a 2 v 2 + a 3 v 3 + a 4 v 4 
or 

Sv s = W -1 Sr + a 4 [v 4 - VW 1 r 4 ]; a 4 arbitrary 

The vector in brackets cannot be zero, since otherwise we obtain the conse- 
quence that s 4 is a linear combination of Si, s 2 , s 3 in contradiction to 
the hypothesis that sx, • • . , s 4 were independent . 1 Since the vector in 
brackets is nonzero, it may be replaced by a unit vector in the same direction 
and, since a 4 is arbitrary, Sv g can now be rewritten as 

Sv s = W -1 Sr + au; a arbitrary (A9) 

where 

u = (v 4 - W _ 1 r 4 )/|v 4 - VW - 1 r 4 | 

In this case there are arbitrarily many satisfactory states corresponding to 
each position deviation, Sr. The union of all such states is then 

S(t) = {x(t) = (Sr,W -1 Sr + au) ; Sr, a arbitrary} (AlO) 

The velocity change required to correct the state (Sr, Sv) to a satisfactory 
state is: 

Av = Sv s - Sv = W 1 Sr + au - Sv; a arbitrary 

Since a is arbitrary in this result, the velocity correction is not uniquely 
specified as a function of the state alone, and a condition that properly spec- 
ifies the value of a must be adopted. In equation (A9) a is arbitrary 
and, therefore, a satisfactory state is achieved independent of the 
u-component of the final velocity. Equivalently, the mission constraints are 
independent of the velocity deviation in the direction u. Consequently, it is 

1 If V 4 - W - 1 r 4 is zero, then the state S 4 can be given as 

34 = (^) = (v) w “ lF4 

That is, s 4 would be a linear combination of s x , s 2 , s 3 with the three 
elements of the vector W _ 1 r 4 as the constants of combination. 
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pointless to make any change in the u -component of velocity deviation. Assum- 
ing that no such change occurs, a must satisfy the condition 2 

_T— 

0 = u Av 


from which 


This gives the result 


a = u T ( -W -1 8r + 5v) 


Av = (I - uir*-)W -1 Br - (i - uu?)Sv 


(All) 


Equation (All) satisfies the definition of a linear impulsive guidance law. 
Its corresponding matrix is then: 


Gi(t) = (I - uu T )W -1 


G 2 (t) = -I + 


— T 

uu 


(A12) 


u = 1 


Case c : n = 5 

The proof of the theorem will not be carried out for this case. The steps 
require only a minor extension of the proof for Case b above. 


Case d: n =6 

This is the trivial case in which every state in N(t) is satisfactory 
and no correction is ever made; that is, 

Gi = G 2 = 0 

2 The condition adopted is equivalent to choosing a such that the magni- 
tude of the velocity correction is minimized. This condition is not required 
by the definition of the guidance law, but is included as a possibility. The 
most general choice of a which satisfies the definition is the form: 

a = (u + p) (-W - Sr + Sv) 

where p is any vector perpendicular' to u. The corresponding guidance law 
matrix will then be 

Gi(t) = -G 2 (t)W _1 
G 2 (t) = -I + u(u + p) T 


This general result points out that the guidance law corresponding to some par- 
ticular four-dimensional set of satisfactory states s(t) is not unique until 
a condition on a is adopted. However, the general result is not of any sus- 
pected practical interest and will not be investigated further. 
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In summary, theorem 1 of the text is now proved. Equations (A 6) an d (A12) 
together with the analogous result for Case c give the expressions for the 
guidance law matrix quoted in the text. 
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APPENDIX B 


UNIQUENESS OF THE GUIDANCE LAW MATRIX 


It is easily shown that the guidance law matrix for the minimum correction 
to a given set of satisfactory states is unique. Let S(t) be a given set of 
satisfactory states and assume that the guidance law exists at t. Let G(t) , 
G*(t) be two guidance law matrices which satisfy the definition in the first 
section and correct any state in N(t) to a state in the given set S(t) . 

Since (l + G*)x(t) is a satisfactory state for any x(t), we have 

G(l + G*) = 0 

G*(l + G*) = 0 


or 

Gj_ = -G2GX* 
Cr ;2 - -G2G2* 
Gj.* = -G 2 *Gi* 


(Bl) 


The set S(t) may be three-, four-, or five -dimensional, and the remainder of 
the proof is carried out according to the dimension of S(t). 


(a) dim [S(t)] = 3 - Theorem 1 gives G 2 , Us* uniquely for this case, 
as 

G 2 = G 2 * = -I 


Then, from the first of equations (Bl) 

G x = Gx* 

that is, G = G*, and the guidance law matrix is unique. 

(b) dim [S(t) ] = 4 . Assuming the condition that the guidance law mini- 
mizes the magnitude of the velocity correction, theorem 1 gives G 2 , G 2 * in 
the form: 

G 2 = -I + uu^ 

G 2 * = -I + u*u* T 

|u| = |u*| = 1 

Inserting these forms in the second equation of (Bl) , obtain 

0 = [u* - (u • u*)u]u* T 
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which, since u, u* are unit vectors, requires that 

u* = ±u 


In either case, 

Cr2 = G 2 * 

and from the first and third equations of (Bl) 

G x = -G 2 G!* = -G 2 *G x * = G x * 

that is, G = G*, and the guidance law matrix is unique. 

( c ) dim [ S ( t ) ] - 5 . If the guidance law is assumed to minimize the 
magnitude of the velocity correction, theorem 1 gives G 2 , G2* in the 
form 


G 2 = -I + U 1 U 1 '*‘ + U2U2^ 

G 2 * = -I + + U2*VL2*^ 

where u x ,U2 and Ui*,U2* are orthonormal pairs of vectors. It is convenient 
to replace these with the equivalent forms: 

G2 = -UsUs^ 

G2* = ~U3*U3* T 


where U-i,Vl Z ,xl 3 and Ui*,u 2 *,U3* are orthonormal triads of vectors. Inserting 
these forms in the second of equations (Bl), obtain 

0 = U3[U3 T - (u 3 . . U 3 *)U3*' 1 '] 

Following the remaining steps of Case b above, obtain 

u 3 * = ±u 3 


G2* = G2 


Gi* = Gi 


that is, G = G* and the guidance law matrix is unique. 



APPENDIX G 


PROPERTIES OP THE MATRIX B 

B = A 4 (t,t 0 ) - B 2 (t,t 0 )Bi 1 (t,t 0 )A 2 (t,t 0 ) (Cl) 

The matrix, Bx(t,t 0 ), has "been assumed nonsingular and the matrices of 
equation (Cl) are defined "by equations ( 6 ) and ( 8 ) of the text. The time 
arguments are dropped in the following. 


B is a Nonsingular Matrix 


Consider the four independent states given hy the columns of 


I o 

Gi(to) ' u 


(G2) 


where u is any nonzero vector, I is the 3 x 3 identity matrix, and Gi(t 0 ) 
is the appropriate submatrix from the guidance law matrix at t 0 . The state 
transition matrix, A(t,t 0 ), is never singular, and, hence, the four columns of 


I Ax A 2 

I A3 A4 



Bx 

A 2 u 

" 

_b 2 

A4U_ 


(C3) 


must also be independent. Suppose, next, that Bu = 0; that is, from equa- 
tion (Cl) : 

A 4 U = B 2 Bx 1 A 2 u 


Note the identity, I = BxBx 1 , from which: 


A 2 u = B xB x A 2 u 


Combining these last two equations yields: 




B -1 AeU 


that is, if 


Bu = 0, 


then the state 



can be given as a linear combination 


of the columns of 



where the constants of combination are the three 


components of the vector Bx 1 A 2 u. However, this contradicts the fact that the 
four columns in equation (C3 ) are independent; consequently, Bu cannot be 
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zero. Since this is true for every nonzero vector, u, it follows that B is 
nonsingular. 

Reduction of B When A 2 is Nonsingular 
The definition of in equations (8) gives 

I = B1B1 = (Ai + AgGi)Bx 
which can he rearranged to give 

I - A2G1BX = AxBi 1 

Assuming A 2 is nonsingular, pre- and post -multiply hy A^ 1 and As, respec- 
tively, and obtain 

I - CriB^As = A 2 1 AiBx 1 A2 (C4) 

Next, rearrange equation (Cl) after introducing the definition of B 2 from 
equations (8) 

B = A 4 [I - G 1 B£ 1 A 2 ] - AsB^As 

The result in equation (C4) is introduced: 

B = [A 4 A2 1 Ai - A 3 ]Bx A 2 (C5) 

T 

The expression in brackets in equation (C5) is [A 2 1 ] , which can be proved 
from the inversion property of the transition matrix 1 (ref. l) 



When the identity I = AA -1 is carried out in terms of its 3X3 submatrices, 
two of the four equations obtained are 

rp T 1 

A 2 Ai - AiA 2 = 0 
A 4 A 1 T - AsA 2 T = I 

From the first of these it follows that A^Ai is a symmetric matrix. Noting 
this fact in a rearrangement of the second equation then yields the indicated 
expression for (A^ 1 )^. Thus, when A 2 is nonsingular, B reduces to 

B = (A ^^ 1 ) T Bi 1 A 2 


1 This inversion property holds for all canonical transformations . 
42 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 
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